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Abstract. We study the shuttling instability in an array of three quantum 
dots the central one of which is movable. We extend the results by Armour 
and MacKinnon on this problem to a broader parameter regime. The results 
obtained by an efficient numerical method are interpreted directly using the 
Wigner distributions. We emphasize that the instability should be viewed as 
a crossover phenomenon rather than a clear-cut transition. 



1. Introduction 

Since the first proposal by Gorelik et al. [T| of the shuttling instability in a 
generic nanoelectromechanical system (NEMS) consisting in a movable single electron 
transistor the phenomenon of the shuttling transport attracted much attention. 
However, until recently the fully quantum theory of the phenomenon was not 
developed. The first quantum mechanical study on a modified setup was accomplished 
by Armour and MacKinnon 2 , closely followed by the work of the present authors 
on the original Gorelik's setup [3] . In this paper we extend the results by Armour and 
MacKinnon. The phase space analysis in terms of the Wigner functions introduced 
in [3] reveals directly the nature of the transport (incoherent tunnelling versus 
shuttling) in different regions in contrast to the indirect evidence used by Armour 
and MacKinnon. Also the new numerical scheme that we use enables us to access a 
wider range of parameters fast and reliably. 

2. Model and method of solution 

Let us consider a simple NEMS consisting of an array of three quantum dots 
(device) connected to two leads 2 . The central dot is assumed to be movable in 
a parabolic potential. The Hamiltonian of the device consists of the mechanical and 
the electronic parts H = iJ mcc h + H c \ where H mc( ± — Uuia'a is the Hamiltonian of 
the harmonic oscillator and H c \ = J2 a /3e{o l c r} \ a ) e aff(x)(P\- We assume strong 
Coulomb blockade regime with no double occupancies so that the vectors \a) with 
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a = 0, L, C, R span the entire electronic Hilbert space of the device. Each matrix 
element e a p{x) is still a full matrix in the oscillator space. e(i) reads explicitly 
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where x = y 2^j(o + <^) i s the position operator, Ay, the device bias, is the difference 
between the energy of the left and the right dot. xq is half the distance between the two 
outer dots and represents the maximum amplitude of the inner dot oscillation. The 
three dots are electrically connected only via a tunnelling mechanism. The tunnelling 
length is given by 1/a and the tunnelling strengths depend on the position operator 
x of the inner grain as tz,(x) — Voe~ a ( Xo+x \ tu(x) — Voe~ a ( x °~ x \ 

The dynamics of the device is described by the generalized master equation 
(GME) formalism The density matrix of the device evolves according to p = 
— i[H, p] + Sp + pd- The first term represents the coherent evolution of the isolated 
device. The coupling to the leads responsible for the electronic transfer from/to the 
device is introduced in the wide band approximation following Gurvitz et al. 0] and 
yields the second term in the equation (in the block notation of 0) 
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with r being the injection rate to the leads. The third term describes the effect of the 
environment on the oscillator, consisting in mechanical damping and random quantum 
and thermal excitation (Langevin force). It reads [2] 



p\i = n(aaJ p — 2a) pa + paa)) — -^(n + l)(a^ap — 2apa^ + pa^a) 
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(3) 



where 7 is the damping rate and n is the mean occupation number of the oscillator 
at temperature T. 

The stationary version of the above GME was solved numerically after the 
oscillator Hilbert space was truncated at sufficiently large N so that all dynamically 
excited states were contained within the basis. Utilizing the decoupling properties of 
the GME in the block notation resulted in the problem of finding a unique null space 
of a (super)matrix with the linear dimension 107V 2 (with up to N = 40). The Arnoldi 
iteration [H] with preconditioning was found to be superior to direct methods, such 
as singular value decomposition or inverse iteration. The preconditioning involving 
the inversion of the Sylvester part of the problem (a fast procedure) is crucial for the 
convergence of the method which is very fast and low memory consuming (the whole 
(super)matrix does not have to be stored in the memory). 



3. Results and discussion 



The resulting stationary density matrix was used to evaluate the mean value of current 
I = r Tr osc pij/j |2] and the phase space distribution of the charged central dot using 
the transformation into Wigner coordinates 0. These quantities as a function of 
the device bias AV and for three different values of the injection rate T are plotted 
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Figure 1. Plot of the current through the triple dot system as a function of the 
device bias for different injection rates T. The other parameters are Vq = 0.5, 
a = 0.2, x = 5, 7 = 0.0125, u = 1. 



in the figures 1 and 2, respectively. It was found in [5] that the triple dot system 
exhibits different regimes of transport at different device biases. The current peaks at 
AV ~ nu (see figure 1) were identified as effects of electromechanical resonances within 
the device. Yet, the different peaks may correspond to different physical mechanisms 
— while the peak around AV w u> is mainly due to the incoherent oscillator-assisted 
tunnelling the peak at AV w 2uj reveals a clear shuttling component. This finding by 
Armour and MacKinnon based on indirect evidence of parametric dependencies of the 
current curves (e.g. the dependence of the current curve on the tunneling length 1/a) 
is confirmed by the direct inspection of the phase space distributions (see the first row 
of figure 2). The half- moon- like shape characteristic for shuttling transport [3] is only 
present around AV^ « 2uj while all other plots show just the fuzzy spot indicative of 
incoherent tunnelling. However, our direct criterion for detecting the shuttling regime 
reveals a close similarity between the resonances. For increasing injection rate we 
can see that the shuttling regime gradually sets in also in the vicinity of the first 
resonance peak. This reveals the crossover character of the onset of the shuttling 
instability found also previously 3 . The sharp transition into the shuttling regime 
reported in semiclassical studies is smeared into the crossover due to the noise present 
in the system and properly accounted for in our approach. 
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Figure 2. Wigner distributions for the central dot in the charged state. Different 
device biases are the points of minimum or maximum current (see figure 1). 
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